Spectral function and quasi-particle damping of interacting bosons in two dimensions 
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agrammatic approaches, and is consistent with the exact Nepomnyashchy identity, which states 
that the anomalous self-energy vanishes at zero frequency and momentum. We recover the correct 
infrared behavior of the propagators and present explicit results for the spectral line-shape, from 
which we extract the quasi-particle dispersion and damping. 
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The excitation spectrum of the weakly interacting Bose 
gas has been studied with field theoretical methods for 
more than half a century [l[ and is qualitatively well un- 
derstood: at high energies the excitations resemble those 
of free bosons whereas at low energies the excitations are 
collective Goldstone modes which result from the symme- 
try broken ground state 0, H, H|. In D = 3, this picture 
has recently been confirmed experimentally using the 
Bragg spectroscopy technique on cold atoms How- 
ever, while perturbative approaches (based on a renor- 
malized T-matrix) to the weakly interacting Bose gas 
have been successfully applied to study the excitation 
spectra and the damping of quasi-particles Q, the self- 
energies in these approaches violate the exact Nepom- 
nyashchy identity [3], which states that the anomalous 
self-energy at zero frequency and momentum vanishes. 
Moreover, the perturbative approach is plagued by in- 
frared (IR) divergences 0, which makes a numer- 
ical evaluation non-trivial already at second order in 
the effective coupling constant. Modern renormaliza- 
tion group (RG) methods [1, [tj [l(| have resolved this 
problem and recovered the correct low-energy structure 
of the self-energies which are in accordance with the 
Goldstone character of the excitations and the Nepom- 
nyashchy identity. The calculations in Refs. [H, 0, US] 
were however limited to the asymptotic regime near zero 
energy and are thus unable to describe both the damping 
of quasi-particles and the crossover from the low energy 
Goldstone-like dispersion to the free particle dispersion 
at higher energies. 

There is a renewed interest in the physics of interacting 
bosons since the transition from weak to strong interac- 
tion can now be studied experimentally via the Feshbach 
resonance technique [ll| . The excitations of the weakly 
interacting gas are qualitatively different from those of 
strongly interacting liquids such as 4 He and the connec- 
tion between the two remains elusive. A consistent strong 
coupling approach, which would allow a description of 
the strongly interacting gas starting from a microscopic 



model is at present not available. In two dimensions, 
the need to go beyond mean field theory is more urgent, 
since the s-wave scattering length vanishes and the usual 
expansion parameter, the gas parameter, is replaced by 
a new emergent parameter }l2j . In this Letter we show 
how the functional RG (FRG) can be employed to calcu- 
late the single-particle spectral density of the interacting 
Bose gas in D dimensions and present numerical results 
for D = 2. Our approach can be extended to include arbi- 
trarily strong local correlations and might offer a possible 
route to the physics of strongly interacting bosons. 

We consider bosons with mass m, chemical potential 
fi and a repulsive contact interaction uq at zero temper- 
ature. The bare Euclidean action is 
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where the spatial integrals should be regularized by 
means of a short-distance cutoff Aq" 1 , which is related 
to the finite extent of the interaction or, for the hard 
core bosons, to the size of particles. Note that the 
model ([T]) depends on two dimensionless parameters 
/to = 2m/LiAo~ 2 and uq — 2mitoA ' )_2 , which give the 
chemical potential and interaction energy in units of 
Ag/2m. The FRG is conveniently formulated in terms 
of the generating functional of irreducible vertex func- 
tions Ta^, <f>], dependin g on an IR cutoff A which even- 
tually will be removed |l3j. The functional T\ is de- 
fined through T A [4>, 4>] = C\[(f>, 4>] + f K 4> K G^ K {K)4> K 
where K = (k,iu) and J K = (2tt)-^ d+1 ^ f d b k doj is 
the integration over momenta and Matsubara frequen- 
cies. Here, C\[4>, <j>] is the Legendre transform of the 
generating functional of connected Green functions and 
Gq \ — iu> — £fe + fi — R\{k) is the inverse free propaga- 
tor with dispersion e& = k 2 /2m. The regulator function 
R\(k) will be specified below. We shall use the following 
non-local potential approximation for T\, 
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Ik 



&Pkuk(K)5p- 



K 



(2) 



A' 



2 



The quantity Spx — Jq <I>q4>k+q — Sk.oPa ^ s ^ ne Fourier- 
transform of px~ Pa where X — (x,t). Here, px — \4>x\ 2 
is the local density and p A is the flowing condensate den- 
sity, which has a finite limit for A — > in the symmetry 
broken state. The non-local character of the potential 
coupling function ua(K) ,which was not included in pre- 
vious FRG approaches 0, [1(3], is essential to describe 
correctly the low energy sector of the theory. Although 
we will in the end calculate the complete momentum and 
frequency dependence of a\(K) and u\(K), we shall also 
employ the following low energy expansion, 

a A (K) « fx(l-X A )+ico(l-Y A ) 

+ e k (Z^ - 1) - (iw)V A , (3a) 
ua(K) w u a (0) = u A , (3b) 

with Xa = p a ua/p- The flow of the parameters Y A , Za, 
Va, Pa> an d u\ will be determined self-consistently. The 
initial values can be read off from the action fl} and are 
given by X Ao = Y A „ = Z Ao = 1, Va = and u Ao = u . 
The initial condition for the condensate density is 
given by the mean field result p Ag = p/uq. The usual Bc- 
liaev Green functions are defined via the matrix equation 
G A \K) = G~\(K) - Sa(A-), where 



Ga(K) 

G q ,a( k ) = 
Va(K) = 
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The normal (E^) and anomalous (E A ) self-energies are 
obtained from an expansion of Eq. ([2]) in 5<fix = 4>x — <fi A 



3, Hj], where we assume a real valued order parameter 
such that Pa — i'A't'A = (^a) 2 - The self-energies are 
obtained from ua{K) and ua(K) via 



E%(K) = aA(K)+p A [uA{0)+UA{K)], 



Hi(K) = p a ua(K). 



(7a) 
(7b) 



Conversely, we may determine (ta{K), which at T = 
has a well defined low-frequency and momentum ex- 
pansion 0], and the non-local coupling function ua(K), 
which will become non-analytic for A — > 0, from a cal- 
culation of the self-energies for all frequencies and mo- 
menta. The effective action ([2]) elegantly divides the 
well behaved parts of the self-energies from the non- 
analytic ones. It further guarantees, since only the flows 
of the [/(l)-invariants ua, p A an d <r\ are retained, that 
the Hugenholtz-Pines relation £^(0) - Ef (0) = /i | 
is obeyed for all A. We use the Litim cutoff function 
R A (k) = (l~S k:0 )(2mZ A )- 1 (A 2 ~k 2 )e(A 2 -k 2 ), 
which has good convergence properties [161 ] and yields 
simple momentum integrals, to derive the flow equations 



for the coupling parameters [^, [nj. We define ka — 
K D [l - r) z A /(D + 2)]/D with K D = 2 1 - D n- D / 2 /T[D/2}, 
where rf A = Ad a In Za is the anomalous dimension of the 
fields, and A a —uap\- We arrive at 
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where T> A (iuj) = Y|w 2 + [e A + Va^ 2 ][£a + Valo 2 + 2A/ 
and £a = (-a/Za- Here we have defined c. 
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The flow of Pa is chosen such that in the expansion 
of Ta in the field 54>x — 4>x — 4>a 110 t erm s linear in 
S(f>x appear 14|, i. e., one always expands around the 
flowing minimum of Ta- For the coefficients entering 



Eq. (|8bj) we have c, 
c («) = 3V A (5? 

(u) 
-2 



(u) 



5i\ 



3e A A A + 6e A A 



4Ai, 



a t 2? A A A + 2A A ) - Y 2 (7A A + lle A ), 
c 2 - = VA[3V A (5tA + AA)~llY 2 ], and t£° = 5V A - 
The flow equations of the derivative terms Ya, Za and 
Va are extracted from the flow equation of 
{K) - (K) = p A u A + a A (K). We find 
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where 

4 V) = V 2 , and 

>) 
c i 



-2Ai, 
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-A 
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)^A, 

a)^a, 



26aA a -^^ a , h 

-r A 2 (e~ A + A A )-e A (? A 
- 2V A {Y 2 + V A (e A + A A )], 4 W) = 3^ A 3 . 
Similar equations were already discussed in Ref . [9j, [l(| 
and we briefly summarize the results in two dimensions. 
At T = the condensate density p A flows to some fi- 
nite limit Pa = o > 0, while the coupling constant ua 
vanishes for A — > which ensures E A _^ (if = 0) = /i 
and T, A ^ (K = 0) = 0, in accordance with both the 
Hugenholtz-Pines relation 0] and the Nepomnyashchy 
identity 0. Furthermore, Ya vanishes for A — * 0, again 
in accordance with exact results Q- Since Va acquires 
a finite value and Za remains finite for A — > 0, the low 
energy modes are phonons with a linear dispersion and 
velocity c = (2mVZ)- 1 / 2 with V = V A =a and Z = Z A =o- 
In Fig. [1] we show flows of the coupling parameters and 
of the dynamical exponent za = 2 — if A — r\ A where 
if A = Ad a In Ya and ?y A = A9a In Za , which conveniently 
illustrates the crossover from the free particle regime with 
z = 2 to the Goldstone regime, where 2 = 1. To work out 
how the crossover appears in the momentum dependence 
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FIG. 1: Typical flows of the parameters Xa, Ya, Za, Va, Pa 
and dynamical exponent za = 2 — r/ A — rf^. The curves are 
calculated for po = 2m/iA^" 2 = 0.4 and uo ~ 2muo = 4. The 
arrows point to the relevant scales. 



of the self-energy, the truncations (|3al3bj) are however 



not sufficient. Moreover, within these truncations one 
obtains G A (K) = liniA->o G A (K) = which violates the 
asymptotic K — > relation [3, [l7 1 



G N (K) ~ — G (if) 



n 2 
mp c 



p (lu 2 + c 2 k 2 ) 



(10) 



where p is the condensate density, p the boson den- 
sity and c the thermodynamic sound velocity. This re- 
lation can be satisfied only within an approach which 
takes account of the non-analytic structure of ua(K) 
Q. We further note that the crossover energy scale 
A c = 47rp°/[TOln( / 9°a 2 )], obtained by Schick [ij] for hard 
core bosons of diameter a, emerges from the z = 2 regime 
of Eqs. (|8bl8a|) in the dilute limit and for large u Aq . 

We now proceed with the evaluation of the self- 
energies. Truncating the momentum dependence of the 
irreducible vertices on the right-hand-side of the flow 
equations in accordance with the truncation of Ta in 
Eqs. ((2]) and (|3al3bp . we find the flow equation of the 
self-energies (the single scale propagators G A ' N (K) are 
defined via G A {K) = -G A (K)d A G-\(K)G A (K) ), 



d A ^(K) = 2u A [ [G N A {Q) + Gi(Q))-Au\pl f [g n a (Q)[G n a (Q + K) + G N A {Q - K) + G N A (-Q + K) 

+2G£(Q - K)] + 2G A (Q) [G A (Q + K) + G A , (Q + K)] } , (11a) 

d A Zi(K) = 2u A [ G N A (Q)~Au\pi [ \g a (Q) [G a (Q + K) + G A (Q - K) + G A (Q + K) + G A (Q - K)] 
JQ JQ L 

+G£(Q) [G% (Q + K) + G N A (Q-K) + 3G£(Q + K)] } . (lib) 

I 



To solve Eqs. (|llalllb[) . we adopt a technique we 
applied previously to the classical </> 4 -model [Tof (see 
la . [l9| for approaches to classical models in their sym- 
metric phase). Instead of attempting to calculate a 
completely self-consistent solution, we use a non-self- 
consistent approach where on the right-hand-side we 
approximate the self-energies entering the Green func- 
tions using Eqs. (|3al3b[) and (|7al7b|) . The full solu- 
tion is expected to appreciably deviate from the non- 
self-consistent solution only for large momenta (at the 
order of A ) or at strong coupling. With Eqs. (|3al3bp 
and the already determined flows of u A , p A ,Y\, Z A and 
Va, we can simply integrate Eqs. (jllalUbp to obtain the 
complete momentum and frequency dependence of the 

self-energies Y, N I A (K) = lim A ^ ^ A M ( k ) ■ Thc result " 
ing expressions for Y, N I A (K) are free from IR divergences 
since the coupling constant u A vanishes in the IR limit. 

An important result of this approach is the correct de- 
scription of the non- analytic structure of u A {K) Q with 



leading terms linear in \lu\ and |fe| (for D = 2). Thc 
Green functions G N (K) and G A {K) are obtained from 
G- l {K) = G^(K) - T,(K) and one can check that 
since p°\u A ^ (K)\ > \u?V + e k /Z\ for small K, G N (K) 
and G (K) have the limiting behavior given in Eq. (jTUJ) 
if we identify Q fP/p =Z and c = {2mVZ)- 1 / 2 . 

The single-particle spectral density A(k,u>) is obtained 
from the imaginary part of the normal propagator via 
A(k,Lu) = —2ImG N (k,iuj — > u + iO) and requires ana- 
lytic continuation to real frequencies. We used the stan- 
dard Pade technique [2(| with 200 points. In Fig. the 
spectral density is plotted as a function of lu > for 
different momenta. One clearly observes a finite peak 
broadening, which grows with increasing momenta. The 
broadening arises from Beliaev damping [||. Because of 
the upward curvature of the dispersion, a momentum- 
and energy-conserving decay of quasi-particles into pairs 
of quasi-particles is allowed. A second order perturbative 
analysis [21( shows that in D dimensions this damping at 
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FIG. 2: Single-particle spectral density A{k,ui) as a function 
of lo, for different values of k, calculated for pLo = 0.15 and 
Mo = 15. The inset shows the quasi-particle dispersion 
which deviates at large k from linearity but is well described 
by the Bogoliubov form E\ = + <?k 2 with renormalized 
velocity c (black dots). The peaks (from left to right) corre- 
spond to k/k c = 0.2, 0.3, 0.4, 0.5 and 0.6, where k c = 2mc. 
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FIG. 3: FRG result for the quasi-particle damping 7 fc. Black 
dots are extracted from the spectral density in Fig. [2] while 
the solid line fits them as 7fc ~ 0.194fc 3 /2mfc c . 

small momenta and for weak coupling has the form 

7 f Kl D (2mp°)^kl- D e D -\ (12) 

where I D = 2~ 4 3 — Kd-i f dx(x~x 2 ) D ^ 1 , fco = 2mco, 
and Co = \J p°uq /m is the mean field velocity. While in 
D = 3 the perturbative expression is independent of cq 
and reproduces the Beliaev result 7^ = 3fc 5 /(6407rm/9°), 
in D — 2 the damping depends explicitly on cq and is thus 
expected to be renormalized. An analysis of the data in 
Fig. [2] reproduces the /c 3 -behavior for small momenta as 
predicted by Eq. (TT2"|) and shown in Fig. 03 However, 
the prefactor a = I D (2mp a y 1 k^ D in Eq.(fT2]l should 
be replaced by a function cv(/io,Mo) of the relevant di- 
mensionless parameters jlo and uq of the model ([1]). For 
example for p. = 0.4 and uq — 4 we obtain c/c ~ 1.01 
and a/a ss 0.915, while the results for p = 0.15 and 
uq = 15 are c/co ~ 0.669 and a/ao « 0.526. The inset 
of Fig. [5] shows the extracted quasi-particle dispersion, 
which is well fitted by E^ = e\+ c 2 fc 2 with renormalized 
velocity c. The deviation of the spectrum from linear- 
ity occurs at the same momentum where the damping 
deviates from the cubic asymptote. 

In conclusion, we have presented FRG results for the 
spectral function of interacting bosons in D = 2 in an 



approach which is consistent with the Nepomnyashchy 
identity T, A (0) = and the Hugenholtz-Pines rela- 
tion £^(0) - £ A (0) = /i. While previous RG calcula- 
tions HI, 0, [HI] were limited to the Goldstone regime, 
our approach allows calculating the entire spectral line- 
shape, including the dispersion and the damping of 
the quasi-particles. Our truncation captures the non- 
analytic structure of the self-energies, described by the 
non-local potential coupling u\(K), which is essential for 
a correct description of the low-energy physics. We are 
currently investigating extensions of our method to in- 
clude arbitrarily strong local correlations. This would 
possibly provide a non-perturbative access to strongly 
interacting bosons. 
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1 



The excitation spectrum of the weakly interacting Bose gas has been studied with field 
theoretical methods for more than half a century [l| and is qualitatively well understood: 
at high energies the excitations resemble those of free bosons whereas at low energies the 
excitations are collective Goldstone modes which result from the symmetry broken ground 
state . In D = 3, this picture has recently been confirmed experimentally using the 

Bragg spectroscopy technique on cold atoms [S]. However, while perturbative approaches 
(based on a renormalized T-matrix) to the weakly interacting Bose gas have been successfully 
app.ied to study the excitation spectra and the damping of auas.partic.es g the self- 
energies in these approaches violate the exact Nepomnyashchy identity [7J, which states 
that the anomalous self-energy at zero frequency and momentum vanishes. Moreover, the 
perturbative approach is plagued by infrared (IR) divergences {2,2], which makes a numerical 
evaluation non-trivial already at second order in the effective coupling constant. Modern 
renormalization group (RG) methods [hJ have resolved this problem and recovered the 
correct low-energy structure of the self-energies which are in accordance with the Goldstone 



character of the excitations and the Nepomnyashchy identity. The calculations in Refs. 
Is, [h3] were however limited to the asymptotic regime near zero energy and are thus unable 
to describe both the damping of quasi-particles and the crossover from the low energy 
Goldstone-like dispersion to the free particle dispersion at higher energies. 

There is a renewed interest in the physics of interacting bosons since the transition from 
weak to strong interaction can now be studied experimentally via the Feshbach resonance 
technique The excitations of the weakly interacting gas are qualitatively different 

from those of strongly interacting liquids such as 4 He and the connection between the two 
remains elusive. A consistent strong coupling approach, which would allow a description of 
the strongly interacting gas starting from a microscopic model is at present not available. In 
two dimensions, the need to go beyond mean field theory is more urgent, since the s-wave 
scattering length vanishes and the usual expansion parameter, the gas parameter, is rep.aced 
by a new emergent parameter [121 ]. In this Letter we show how the functional RG (FRG) 
can be employed to calculate the single-particle spectral density of the interacting Bose gas 
in D dimensions and present numerical results for D = 2. Our approach can be extended to 
include arbitrarily strong local correlations and might offer a possible route to the physics 
of strongly interacting bosons. 

We consider bosons with mass m, chemical potential \i and a repulsive contact interaction 
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Mo at zero temperature. The bare Euclidean action is 



S[if>,iJ)] = / dPxdr 



__2 

^~2^ -^ + y(#) 2 



(1) 



where the spatial integrals should be regularized by means of a short-distance cutoff A , 
which is related to the finite extent of the interaction or, for the hard core bosons, to 
the size of particles. Note that the model (pQ) depends on two dimensionless parameters 
Jxq = 2m/iAo~ 2 and uq = 2mwoA :? ~ 2 , which give the chemical potential and interaction energy 
in units of !\\j1m. The FRG is conveniently formulated in terms of the generating functional 
of irreducible vertex functions T A [(j), (/>], depending on an IR cutoff A which eventually will be 
removed [3]. The functional T A is defined through 0] = £ A [(j),<f)] + J K 4>kGq \{K)4>k 
where K = (k,icu) and J K = (2n)~ ( ' D+1 ^ J d D k du is the integration over momenta and 
Matsubara frequencies. Here, C A [<f), <P] is the Legendre transform of the generating functional 
of connected Green functions and Gq \ = iu — + p, — R A (k) is the inverse free propagator 
with dispersion e k = k 2 /2m. The regulator function R^{k) will be specified below. We shall 
use the following non-local potential approximation for Ta, 

T A = / 4>k(Tk{K)(I)k + \ [ 6p K u A (K)5p„ K . (2) 

JK 1 JK 

The quantity 5px = Jq 4>q<Pk+q — ^k,oPa is the Fourier-transform of px — p\ where X = 
(x, r). Here, px = \4>x\ 2 is the local density and p A is the flowing condensate density, which 
has a finite limit for A — > in the symmetry broken state. The non-local character of the 

B)tential coupling function u\(K), which was not included in previous FRG approaches 
is essential to describe correctly the low energy sector of the theory. Although we 
will in the end calculate the complete momentum and frequency dependence of o~\(K) and 
u\(K), we shall also employ the following low energy expansion, 

a A (K) w fi(l-X A )+iw(l-Y A ) 

+ e k (Z A l -l)-(tiu) 2 V A , (3a) 
u A {K) w u A (0) = m a , (3b) 

with X A = p A u A / fi. The flow of the parameters Y A , Z A , V A , p A , and u A will be determined 
self-consistently. The initial values can be read off from the action (CQ) and are given by X Ao = 
Y Ao = Z Ao = 1, V Ao = and u Ao = u . The initial condition for the condensate density 



3 



is given by the mean field result p° = fi/uo- The usual Beliaev Green functions are defined 
via the matrix equation G^{K) = Gq A (K) — Ha(K), where 

, G N JK) Gi(K) , 
G 0)A (-K) 

S A if = M ' M ; . (6) 



The normal (E^j and anomalous self-energies are obtained from an expansion of 

nn 

Eq. (|2J) in <50x = 4>x — 4>\ L3J| , where we assume a real valued order parameter such that 
Pa = 0a0a = (0a) 2 - The self-energies are obtained from cr\(K) and u\(K) via 



(#) = <x A (if) + p A K(0) + u A (2iQ] , (7a) 
E£(K) = pW(K) • (7b) 

Conversely, we may determine a\(K), which at T = has a well defined low- frequency 
and momentum expansion [7], and the non-local coupling function ua(K), which will 
become non-analytic for A — > 0, from a calculation of the self-energies for all frequen- 
cies and momenta. The effective action ([2]) elegantly divides the well behaved parts 
of the self-energies from the non-analytic ones. It further guarantees, since only the 
flows of the [/(l)-invariants «*, p A and a\ are retained, that the Hugenholtz-Pines re- 
lation E£(0) - ££(0) = p m is obeyed for all A. We use the Litim cutoff function 
RAk) = (1 — 5fc,o) (2^^a) _1 ( A 2 ~~ ^ 2 ) ® ~ ^ 2 )? which has good convergence properties 



16l | and yie 



ds simple momentum integrals, to derive the flow equations for the coupling 
parameters ^Q. We define k a = K D [1 - r] z J(D + 2)]/D with K D = 2 1 - d tj- d ' 2 /T[D /2], 
where rj A = A<9 A lnZ A is the anomalous dimension of the fields, and Aa = u\p°\- We arrive 
at Q 

A D+Y „ /■ J , 3 „(P), ,2n 
n=0 iv 

2mZ A / 2^^pf7i^) 



2 A^a / da; \- 

2mZ A J 2vr ^ D A za; 
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where V A {iui) = Y A u 2 + [e A + V A u 2 ][e A + V A u! 2 + 2A A ] , and e A = e A /Z A . Here we have defined 



» 



A 2 M 
~ A A' c l 



FA(2eA + A; 



■ Y A , and c 2 



to 



V 2 . The flow of p A is chosen such 



that in the expansion of r A in the field 5<px — <Px ~ 0a no terms linear in 6<px appear [14| , 
i. e., one always expands around the flowing minimum of r A . For the coefficients entering 
Eq. §b|) we have c Q u) = 5e A + 3e A A A + 6e A A 2 + 4A|, cf ] = 3V A (5e| + 2e A A A + 2A 2 ) - 



F A 2 (7A A + lle A ) 



» 



V A [3V a (5e A + A 7 



11F A ], and c 3 



(u) 



5V A . The flow equations 



of the derivative terms Y A , Z A and V A are extracted from the flow equation of 
Ef(K) - Ei(K) = p° A u A + a A (K). We find 



A D+1 i 



0„,2t/ 1 K A 



2m Z * 



271 h> v ^ iuj ) 



d A Z A 



d A V A 



4/M 



A D+1 Kj 
1 2mD 



do; 



2%V 2 A (iuj) 



2 iV fi, A 



2mZ 



A 



duj 
2^ 



n=0 



where c, 



to 



;2 



2e A A A - 2A 

» 



,(2/) 



A) u l 



r A 2 + 2(e~ A -A A )^ A , 4 J 



to 



V^, and 



(9a) 
(9b) 
(9c) 
-YHe A + 



A A ) - e A (e A + 2A A )V A , cf> = 2V A [Y 2 + V A (~e A + A A )], a 



(v) 



3Vl 
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FIG. 1: Typical flows of the parameters X A , Y A , Z A , V A , p\ and dynamical exponent z A = 
2 — r] z A — t/ a . The curves are calculated for fi$ = 2?tt,//Aq 2 = 0.4 and uo = 2muQ = 4. The arrows 
point to the relevant scales. 

Similar equations were already discussed in Ref. 0, [h]] and we briefly summarize the 
results in two dimensions. At T = the condensate density p° A flows to some finite limit 
p A=0 > 0, while the coupling constant u A vanishes for A — > which ensures T, A ^ (K = 
0] = p and T, A ^ (K = 0) = 0, in accordance with both the Hugenholtz-Pines relation 
4j and the Nepomnyashchy identity [7]. Furthermore, Y A vanishes for A — ► 0, again in 
accordance with exact results [7|. Since V A acquires a finite value and Z A remains finite 
for A — > 0, the low energy modes are phonons with a linear dispersion and velocity c = 



^raV^Z*) -1 / 2 with V* = V\ =0 and = Z\ =0 . In Fig. [I] we show flows of the coupling 
parameters and of the dynamical exponent z k = 2 — rj\ — n]\ where rf k = A<9a In Y\ and 
r] A = Ac^ln^A, which conveniently illustrates the crossover from the free particle regime 
with z = 2 to the Goldstone regime, where z — 1. To work out how the crossover appears 
in the momentum dependence of the self-energy, the truncations (I3af3bp are however not 
sufficient. Moreover, within these truncations one obtains G A (K) = liniA-^o G A (K) = 
which violates the asymptotic K — > relation 

G"(K) ~ -G\K) ~ mp ° c2 5 1 , (10) 

p (ur + crAr) 

where p° is the condensate density, p the boson density and c the thermodynamic sound 
velocity. This relation can be satisfied only within an approach which takes account of 
the non- analytic structure of u\(K) p|. We further note that the crossover energy scale 
A c = ATTp°/[m ln(p°a 2 )], obtained by Schick [l2] for hard core bosons of diameter a, emerges 
from the z = 2 regime of Eqs. (I8b|8al) in the dilute limit and for large u\ . 

We now proceed with the evaluation of the self-energies. Truncating the momentum de- 
pendence of the irreducible vertices on the right-hand-side of the flow equations in accordance 
with the truncation of T\ in Eqs. (j2J) and ( I3all3b[) . we find the flow equation of the self-energies 
(the single scale propagators G A ' N (K) are defined via G^{K) = —G\(K)d\GQ\(K)G\(K) 

), 



d A ^(K) = 2uJ {G%(Q) + Gi(Q)}-4ulpl j {gZ(Q)[G%(Q + K) + G N A {Q ~ K) + G%(-Q + K) 

J Q J Q 

+2Gi(Q - K)] + 2Gi(Q) [G A K {Q + K) + G%(Q + K 



11a 



<9 a £aW = 2« A / G N K {Q)-Au\ P l [ {g»(Q)[G%(Q + K) + G N k {Q - K) + G A k {Q + K) + G A (Q - K) 



+G£{Q) [G% {Q + K) + G% {Q-K) + 3Gi{Q + K 



To solve Eqs . (fTTalll 1 bl) . we adopt a technique we applied previously to the classical 



model |15fl (see 



18 



19( | for approaches to classical models in their symmetric phase). Instead 



of attempting to calculate a completely self-consistent solution, we use a non-self-consistent 
approach where on the right-hand-side we approximate the self-energies entering the Green 
functions using Eqs. (I3a|3bp and (I7a|l7bp . The full solution is expected to appreciably deviate 
from the non-self-consistent solution only for large momenta (at the order of A ) or at strong 



'lib 



coupling. With Eqs. fl3all3bl) and the already determined flows of u\, p\, Ya, and V\, 
we can simply integrate Eqs. flllafllbl) to obtain the complete momentum and frequency 
dependence of the self-energies %n/A(k) = lim A _ >0 S^ M (K). The resulting expressions for 
Y> N I A (K) are free from IR divergences since the coupling constant ua vanishes in the IR 
limit. 

An important result of this approach is the correct description of the non-analytic struc- 
ture of u\(K) [tJ with leading terms linear in \u\ and \k\ (for D = 2). The Green functions 
G N (K) and G A (K) are obtained from G' 1 ^) = G \K) - E(iT) Q and one can check 
that since p°\u A _ (K)\ > \uj 2 V* + e fc /Z*| for small K, G N (K) and G A (K) have the limiting 
behavior given in Eq. (TTTJT) if we identify p°/p = Z* and c = (2mV^Z^ 1 ^ 2 . 
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FIG. 2: Single-particle spectral density A(k,uj) as a function of u, for different values of k, calcu- 
lated for po = 0.15 and uq = 15. The inset shows the quasi-particle dispersion which deviates at 
large k from linearity but is well described by the Bogoliubov form E 2 = e\ + c 2 k 2 with renormal- 
ized velocity c (black dots). The peaks (from left to right) correspond to k/k c = 0.2, 0.3,0.4,0.5 
and 0.6, where k c = 2mc. 
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FIG. 3: FRG result for the quasi-particle damping 7^. Black dots are extracted from the spectral 
density in Fig. [2j while the solid line fits them as 7^ ~ 0.194/c 3 /2mfe c . 

The single-particle spectral density A(k,u) is obtained from the imaginary part of the 
normal propagator via A(k,u) = —2lmG N (k,iuj — > uj + iO) and requires analytic contin- 
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uation to real frequencies. We used the standard Pade technique [2Cj with 200 points. In 
Fig. [2] the spectral density is plotted as a function of uj > for different momenta. One 
clearly observes a finite peak broadening, which grows with increasing momenta. The broad- 
ening arises from Beliaev damping [6] . Because of the upward curvature of the dispersion, 

a momentum- and energy- conserving decay of quasi-particles into pairs of quasi-particles is 

I I 

allowed. A second order perturbative analysis [2l| shows that in D dimensions this damping 
at small momenta and for weak coupling has the form 

7 f « I D {2mp°)^kl- D k™-\ (12) 

where I D = 2~ 4 3^~ Kd-i Jq dx(x — x 2 ) D ' 1 , k = 2mco, and Co = a/ p°u /m is the mean field 
velocity. While in D = 3 the perturbative expression is independent of Co and reproduces 
the Beliaev result 7^ = 3/c 5 /(6407rmp°), in D = 2 the damping depends explicitly on Co 
and is thus expected to be renormalized. An analysis of the data in Fig. [2] reproduces the 
fc 3 -behavior for small momenta as predicted by Eq. ( TT21) and shown in Fig. [3j However, the 
prefactor a = I D (2mp°)~ 1 kl~ D in Eq. ffl2|) should be replaced by a function a(p , u ) of the 
relevant dimensionless parameters jlo and Uq of the model (Q. For example for p Q = 0.4 and 
uq = 4 we obtain c/cq ~ 1.01 and a/a^ ~ 0.915, while the results for ft = 0.15 and u = 15 
are c/co ~ 0.669 and a/ao ~ 0.526. The inset of Fig. [2] shows the extracted quasi-particle 
dispersion, which is well fitted by E\ = t\ + c 2 k 2 with renormalized velocity c. The deviation 
of the spectrum from linearity occurs at the same momentum where the damping deviates 
from the cubic asymptote. 

In conclusion, we have presented FRG results for the spectral function of interacting 
bosons in D = 2 in an approach which is consistent with the Nepomnyashchy identity 
T, A (0) = and the Hugenholtz-Pines relation S^fO) — S A (0) = u. While previous RG cal- 

nnn 

culations [8|, [9J, HDJ were limited to the Goldstone regime, our approach allows calculating the 
entire spectral line-shape, including the dispersion and the damping of the quasi-particles. 
Our truncation captures the non-analytic structure of the self-energies, described by the non- 
local potential coupling u\(K), which is essential for a correct description of the low-energy 
physics. We are currently investigating extensions of our method to include arbitrarily 
strong local correlations. This would possibly provide a non-perturbative access to strongly 
interacting bosons. 

We thank A. L. Chernyshev for detailed discussions on Beliaev damping, H. O. Jeschke 
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